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Abstract: 
A recursive method is derived to calculate all eigenvalue correlation functions of a 
^ ■ random hermitian matrix in the large size limit, and after smoothing of the short scale 

^ I oscillations. The property that the two-point function is universal, is recovered and the 

three and four-point functions are given explicitly. One observes that higher order corre- 
lation functions are linear combinations of universal functions with coefficients depending 
on an increasing number of parameters of the matrix distribution. 
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In a recent article [Q, Brezin and Zee have calculated explicitly correlation functions 
of eigenvalues of a class of stochastic hermitian matrices of large size A^. They have found 
that some statistical properties of the eigenvalues are universal in the large A^ limit, and 
can thus be obtained from the correlation functions of the gaussian model. More precisely, 
they have discovered that the two-point correlation function, after smoothing of the short 
scale oscillations, is universal while all other correlations vanish at the same order. 

This property can be related to a renormalization group analysis [|],^ which has shown 
that the gaussian model is a stable fixed point in the large N limit. 

Their analysis is based on the, by now standard, method of orthogonal polynomials. 
An essential ingredient in the final answer is a proposed ansatz for an asymptotic form 
of the orthogonal polynomials Pn in the limit N -^ oo and N — n finite. In ref. [Q the 
ansatz is verified in the case of even integrands, and only up to an unknown function. 

Here, we propose a direct proof of the ansatz, using a saddle point method, which does 
not depend on the parity of the integrand, and which allows to determine the previously 
unknown function. Moreover, using a completely different approach ||^], we present a 
recursive method to evaluate all smoothed correlation functions at leading order. We give 
the three and four-point functions explicitly. 

1. Correlation functions of eigenvalues 

Let us first explain the problem and recall the method used in ref. |^] to explicitly 
evaluate the eigenvalue correlation functions. 

We consider N x N hermitian matrices M with a probability distribution of the form: 

P(M) = -e ^ (1.1) 

where V{M^ is a polynomial, and Z the normalization (i.e. the partition function). We 
want to derive the asymptotic form for A^ large of various eigenvalue correlation functions. 
All can be derived from the correlation functions of the operator 0{X): 



1 1 ^ 

O(A) = -tr5(A-M) = -5^5(A-^0 (1.2) 



i=l 



(the /ii being the eigenvalues of M). Indeed, for any set of functions (/i, ..., /fc), one has: 
^ (tr/i(M) . . .tr/fc(M)) = j dX,h{X,) . ..dXkfk{Xk) (0(Ai) . . .0(A,)) . 



The correlation functions of the operator 0{X) can in turn be expressed in terms of the 
partiaUy integrated eigenvalue distributions like p(A) the density of eigenvalues, which is 
the probability that A belongs to the spectrum of M, P2(A, /u) the probability that A and fx 
are simultaneously eigenvalues of M and more generally Pn(Ai, . . . , A^) the probability that 
Ai ...A„ are simultaneously eigenvalues of M. For Ai 7^ A2 . . . 7^ A„ (else some additional 
contact terms have to be added) we find 

{0{X^)0{X2) . . . 0{K)) = A._^p^(Ai, . . . , A.). 

Actually the interesting functions are not directly the 0{X) correlation functions, but 
their connected parts. Indeed, at leading order, when A^ -^ 00, we have the factorization 
property 

(0(Ai)0(A2)...0(AJ)~ n <^(^^))' 

1=1, ...,n 

and thus no new information can be obtained from the complete n-point function. The 
connected function which will be denoted 

(O(Ai)...O(A.)),„„„ = O.(Ai,...,A0 

is only of order 1/A", and thus is a subleading contribution. The method of orthogonal 
polynomials allows to determine directly all these connected functions from only one aux- 
iliary function k(A, fx) = K{p, A) (see refs. [|l],|^ or appendix 2 for details). For instance we 

have: 

Oi(A) = p(A) = /.(A,A), 

02(A,^) = -«:2^A,^), 
03(Ai, A2, A3) = 2a^(Ai, A2)«:(A2, A3)/€(A3, Ai), 
C4(Ai, A2, A3, A4) = -2k(Ai, A2)k(A2, A3)k(A3, A4)k(A4, Ai) 

- 2k(Ai, A3)k(A3, A2)k(A2, A4)k(A4, Ai) 

- 2k(Ai, A2)k(A2, A4)k(A4, A3)k(A3, Ai), 
and analogous expressions for larger values of n: 

1 " 

permutations a i=l 



(each product of k appears exactly n times in the sum and this counting factor is canceUed 
by the 1/n factor). O^ may also be written: 

n 
cyclic permutations a i=l 

The function k is simply related to the polynomials Pn orthogonal with respect to the 

-NV{X) 
measure dAe : 

k(A,^) oc — ^— ^- exp[-{N/2){V{X) + V{fi))]. (1.3) 

The asymptotic evaluation of correlation functions is reduced to an evaluation of the func- 
tion k(A, fi) and thus of the orthogonal polynomials Pn(A) in the peculiar limit: A^ large, 
N — n finite, and A G [a, 6] ([a, b] being the support of p(A)). 

The ansatz proposed in ref. |^ in the case of even potentials V{M) (for which b = —a) 
was: 

P.(A) (X e^^^^^/^ —1= cos (iVC(A) + (iV - nMX) + x(A)), 

where 

A =a cos (p , 

/(A) = asirnp, 

Ac(A) = -np{X), 

x(A) remaining undetermined, except in the gaussian and quartic cases, for which: x(A) = 
ip/2 - n/4. 

We will prove below that this ansatz is still true for any V, and that x is always of 
the form x = V^/2 + const. The method which leads to the proof is also interesting in itself 
because it uses some general tools of the saddle point calculations of the one-matrix model 

m- 



2. Orthogonal Polynomials 

Let us consider the set of polynomials Pn {n is the degree), orthogonal with respect 
to the following scalar product: 

f -NV(X) 
{Pn ■Pm) = Snn.= I dXe P^(A)P^(A). (2.1) 



Remarkably enough an explicit expression of these orthogonal polynomials in terms of a 
hermitian matrix integral can be derived (see appendix 1): 

. ^ -1/2 /■ 2 -mrV(M) ^ ^ , ^ 

Pn{X) = Zjj^ / d" Me ^ ^ det(A - M), (2.2) 



where M is here a n x n hermitian matrix, and the normalization Z^+i the partition 
function corresponding to the n + 1 x n + 1 matrix integral. We will use this expression 
to evaluate Pn in the relevant limit, i.e A^ ^ 1 and N — n = 0{1) by the steepest descent 
method. 

Let us first introduce the matrix integral: 



where 

V{z) = V{z)-hln{X-z). 

With these definitions P„, is proportional to Z{g = n/N, h = 1/N, A), and thus we need Z 

or equivalently the free energy F^ = — n~^ \nZ for h and g — 1 small. 

As we are interested only in the A dependance of F, let us differentiate F with respect 

to A: 

—- = uj{z = \g,h,X) 

dX g 

where uj{z) is the resolvent: 

^(^) = -(tr^— y (2.3) 

n \ z — M / 

Since we want the asymptotic expression of uj{X)'n?h/ g, we need u up to the order 
1/n. It is known from the random matrix theory [^, that F^, and also uj have an expansion 
in powers of l/n^ which is the topological expansion. At order 1/n, only the contribution 
of the sphere is required. We thus replace u by its dominant contribution obtained by the 
saddle point method. 

With this approximation, w may be written: 

-(-) = - E ^ 

n ^-^ z — Xi 

1=1,...,?! 

where the Xi are the eigenvalues verifying the saddle point equation: 



g w^. Xi-Xj 



i.e. they extremize the integrand of 

Z{g, h,X)^ j\[ e~^^^^') (A, \[{X, - X 



In the large n Umit the Xi are distributed along an interval [a, h] with a continuous density 
p(A). Then: 

uj{z)= f di^pii^)^— (2.4) 

J a Z~ ^i 

and the saddle point equation becomes: 

-V'{p) =uj{fi + iO)+uj{fi-iO) for // G [a, 6] (2.5) 

Note an important property of this equation: At a and b fixed it is linear and therefore 
the derivatives of uj with respect to g or h will also satisfy a linear equation. 

At leading order we introduce the resolvent u!q{z) = u!{z, g = l,h = 0), and write: 



gu{z, g, h, A) = uo{z) + {g - l)^g{z, A) + hflh{z, A) + 0(1/ 

where we have defined the two functions: 

dguj 



n 



' dg 



O 



h 



9=0,h=l 9h 



g=0,h=l 



The function uJo{z) is the resolvent of the usual one-matrix model, and from eqs. (^ 
and (^]^) we obtain: 

uJo{z ± zO) = -V'{z) =F mp{z) for z G [a, h]. 

As we noted above, the two functions O^ and O^ obey linear equations, obtained 
by differentiation of the equation ( p.5| ) satisfied by w. Following a method introduced 
in ref. Q we can then easily determine them from their analyticity properties, and the 
boundary conditions. 0^(2;) verifies the linear equation: 

Og(z + iO) +Og(2-zO) =0 
5 



and behaves as 1/ z when 2 ^ oo, and as 1/ \/z — a\/z — b near the cut end-points a and 
6, because u behaves as \/z — a^ z — b. These conditions determine Qg uniquely: 

y/{z-a){z-b) 
The same method apphes to flh which satisfies 

Qhiz + iO) + Qhiz - iO) = 1/(A - z), 

and behaves hke O^ ~ 0(l/z^) for z large, Qh ~ I/V-2 — flV-z — ^ near a, 6, and is regular 
near z = X. It follows: 

Then if we set 2 = A: 

We now have the necessary ingredients to determine dF/dX: 

dF 
-n^-- = Nuo{z = A) + (n - A^)03(2 = A) + 0/,(2 = A) + 0(1/Ar). 
C/A 

We still have to integrate all these terms with respect to A. In order to integrate 

"^0 = V /2 — iivp, we introduce a function ({X) 

-A 



C(A) = -TV / dAV(A'). 



We also need the integral of [(A — a){b — A)]~^/^. For this purpose, we parametrize 
A = I (a + 6) — ^(6 — a) cosi/?, so that the integral of [^{b — a) siuLpj is simply </?. Finally, 
the result takes the form: 

-n^F{X ± zO) = lNV{X) ± iNC{X) T i{N - n)ip ± i^ip 
— ^ In a/(A — a) (6 — A) + const . 

Since P„ is a polynomial, we have Pn(A) = ^ [-Pn(A + zO) + Pn{X — zO)], and therefore: 



where /(A) = a/(A — a){b — A). The constant factor ^/2/t: is fixed by the condition that 
k(A, A) = p(A). In the case of even potentials V, parity considerations imply that the 
arbitrary constant phase is const = — 7r/4. Indeed, we have ({a) = 0, ({b) = — tt, f{a) = 0, 
Lp{b) = TT, ( we have /(a) = f{b) = 0, but for e a small positive number ■.f{a+e)/f{b—e) = 1) 
and thus: 

Pn{b-e)/Pn{a + e) ~ (-l)'"+i tan (const). 

For general potentials the constant phase remains undetermined at this order, but we note 
that the general form of Pn does not depend on the parity properties of V. 

3. Connected Correlation Functions 

From this asymptotic expansion of Pn one can now derive the kernel k(A,/u), and 
then the connected correlation functions, in the large A^ limit. The authors of ref. [0 have 
calculated some correlation functions in two regimes: short range correlations (A^ — Aj) ~ 
1/A^, and mesoscopic correlations (Ai — Aj) ^ 1/A^. Note that the polynomials Pn and 
Pn-1 oscillate at a frequency of order A^ (which corresponds to the discrete spectrum of 
a matrix of size A^ finite), and therefore, all the correlation functions will present such 
oscillations. 

In the short distance regime, these oscillations give the dominant behaviour, and 
eq. (p.6|) leads to: 



sm[2nN{X-^)p{X)] 
"^^'^^ 2.Ar(A - ^) 

and all connected correlation functions follow. Note that k being of order 1/A^, the con- 
nected microscopic n-point function will be of order 1/N^. All this is studied in detail in 
ref. |jl|], and we will now concentrate our attention on the mesoscopic case. 

In the regime A^ — Aj ^ 1/-^, it is interesting to consider smoothed functions, defined 
by averaging the fast oscillations. For instance we find that: 

\n (\ M -1 1 l-cos(^cosV; 

[U2 [A, /ij smooth — ^,.9 2 7"^ ^2 '■ '■ 1 ' ^'^- 

^A/^TT^ (A — ^)^ sm(/?sm'i/; 

where 

^ a + b b — a a + b b — a , ,^^, 

A = ^ ^cos(^, ^=^ —cosip. (3.2) 

Brezin and Zee noted (ref. [jl|) that the smoothed higher order n-points correlation 
functions vanished identically at the order l/N"^ for n > 2. Indeed, we will prove below 
that they are of order 1/A^^*^"~^-*, and give a recursive method to compute them with the 
help of loop-correlators. 



4. The n-Loop Correlation Functions 

Let us introduce the functions: 

.„(........ '„) = iV"-(tr^><...xtr-l^)^^____ (4.1) 

They are related to the previous correlation functions by the relations: 

a;„(zi, . . . , zn) = iV2"-2 f rr -Ah^ 0^(Ai, . . . , A„) (4.2) 

J -^ Zi — \i 
1=1 

ujji is called the n-loop correlation function, because it is the Laplace-transform of the 
partition function of a discrete random surface limited by n loops (the Zi are conjugated 
to the lengths of the loops). This remark allows to understand the topological origin of 
the factor A^^"^~^:indeed, 

the Laplace-transform of the complete n-point correlation function which is of or- 
der 1, would be the partition function of every surface (not necessary connected) with n 
boundaries. Each surface contributes with a topological weight A^-^ where x is its Euler 
character. The leading term is the most disconnected one, with x = n (indeed, such a 
surface is made of n discs, each of them having x = 1), while the connected term has 
X = 2 — n (it is a sphere (x = 2) from which n discs have been removed). Therefore, the 
relative contribution of the connected part to the complete n-loop function is A^2-2n 

Remark that relation ( [4.2| ) can be inverted: u^ is analytical except when some of the 
Zi belong to the interval [a, 6]. On can then be expressed in terms of the differences of Un 
between opposite sides of the cut. For instance: 

p(A) = ^{u{X + zO) - u{\ - zO)) 

ZZTT 

^2(A, ^i) = 777-^^77^ (^2 (A + iO,ii + iO) - uj2{X + iO, ^ - zO) - uj2{X - zO, ^ + zO) 

[ZmJM )^ 

+ U2iX-iO,n-iO)). 
and for general n: 

^n{Xi. . . . An) = j;^ (^) J2 (-l)^''+-+'"^cc;^(Ai + eizO,...,A^-fe„zO) (4.3) 

Note that these functions are directly the smoothed correlation functions, since we first 
compute uJn in the large N limit at complex arguments (which suppresses the oscillations), 
and then take the discontinuities along the cut. 



In order to compute cOn, we consider the partition function: 

^2p r -mrV(M) 

Z = e^ ^ = / dMe ^ ' 



with V{z) = Yl'kLidkZ^ /^^ ^^d define the loop insertion operator: 



^ Y^l_ (4 4) 



5Viz) ^ ^^=+1 dgk 

Note that with this definition 

SV{z') _ 1 _ 1 

dV{z) z — z' z 

The LJn are obtained from the free energy F by the repeated action of this operator: 



(4.5) 



Unizu...,z^) = —— —■■■—— -—F, (4.6a) 

dV{zi) dV{Zn) 

W^_l(2:i,...,2:„_l). (4.66) 



6V{z 

It is not necessary to calculate the free energy F, since we already know the one-loop 
function u!{z). We have already emphasized that u!{z) satisfies a linear equation (eq. ( |2.5D ), 
and thus, all its derivatives satisfy the same linear equation, with a different l.h.s. . Then, 
analyticity properties, and boundary conditions determine the form of cOn- 
The linear equation for u is: 

uj{X + iO)+uj{X-iO) = V'{X). (4.7) 

By a repeated action of the loop insertion operator, we obtain: 

uj2{X + iO,z)+uj2{X-iO,z) = -- — , (4.8) 

{z- xy 

and for n > 2: 

Unix + iO, Z2,.. ., Zn) + UJniX - iO, Z2, ■ ■ ■ , Zn) = . (4.9) 

The function uj{z) has the form 

u;{z) = \ (V\z) - M{z)^{z-a){z-h)^ 



where M{z) is a polynomial such that uj{z) ^ 1/z when z ^ oo, and therefore: 

M(.) ■ ^''^' 



y/{z-a){z~b)J ^ 

Then 

p(A) = ^M(A)v/(A-a)(6-A) = ^M(A)^sin(^. 



where we have used the parametrization (|3.2| ). 

The two-loop function is also completely determined by the linear equation ( [4 .81 ) and 
boundary conditions: 

11 / {x-yf -{x-a){x-b)-{y-a){y-b) 
^{x-yY\ ^^^^a){x-b)^{y-a){y-b) 

and thus, in agreement with eq. (|3.1|) : 

n (\ )= ^ ^ 1 - cos y? cos V> 

i^2lA,^j 2Ar%2(^_^)2 sin(/psini/; ' 

The other loop functions all satisfy an homogeneous equation, and can be written: 

/ n \ 2n-3 

w„(Ai, . . . , A„) = TT ^ Vr„(Ai, . . . , A„), 

with now A„ = ^ ^ cos ifn , and where the Wn are some symmetric polynomials 

of degree less than 2n — 5 in each A^, which are no longer determined by the boundary 
conditions. It is necessary to directly use the recursion relation ( [1.66| ). Since UJ2 depends 
on the potential V{M) only through a and 6, we need the actions of loop-insertion operator 
on a and 6, for instance da/dV{z), d'^a/dV{z) SV{z').... For this purpose, we introduce the 
following moments of the potential: 

Mu = (b dz ^ ^ 



2iiT J (z-a)'' ^(z-a){z-b) 
Jk = -77— <? dz 



2i7T J {z- bf ^^z-a){z-b) 

(the integration path turns clockwise around the cut [a, 6]). 
They are such that: 

M{z) =J2Mk+i{z -a)'' = Y, Jk+i{z - by 



10 



The Mk and Jk are linearly related to the coefficients of the potential V , and if y is a 
polynomial of degree v, M is of degree v — 2, and there are only v — 1 independent coefficient 
among the Mk and Jk- Note also, that if V is even, we have Mfc = —[—l)^Jk- 

The cut end-points a and h depend on V through the conditions that (see ref. ): 

1 /.. ^'w 



2m J ^(z-a){z-b) 
1 I,, ^V'(.) 



^in J ^^z-a){z-b) 



It follows (using ( [4 .51) and performing the contour integrals): 

5a 1 1 1 



Sb 1 1 1 



5y(z) Ji 2 - 6 ^(z-a){z-b) 

In order to determine the higher order derivatives of a and b, we need to differentiate the 
coefficients Mk and Jk- 



5Mk 2k + 1 5a 



K-ir^) 



5V{z) 2 5V{z) V 

1 ^6 / Ji J^_ _ '^ Mi+i 

2 5V{z)y{b-a)>' {z-a)>' ^^ {b - a^-^ 

and analogous formulae are obtained for the Jk by the exchange a ^->- b. 

With these tools, we can now determine the uJn recursively. Let us for instance write 

W3: 

, a — b 

uj2,[x,y,z) 



a/(x - a)(x - 6)A/(y - a){y - b)^{z - a){z - b)j 

' 1 1 

X ( -^{-^ - ^){y -b){z-b)-—{x- a){y - a){z - a] 



For c<;4 we first define the polynomials 

Q{xi,a) = Y[ ( 



^Xi — a] . 

i=lA 



11 



With this notation 

uJA{x^) = — — — —^ ( -3(6 - a)— |Q(x„ h) - 3(a - h)^Q{x,, a) 



1 

1 

-'1 



+ 3772'5(a;i,^) 



+ 3-p2(5(xi,a) 



<.-)(e:.^)- 



+ [(xi - a)(x2 - a)(x3 - 6)(x4 - 6) + 5 terms] 

where the last additional terms symmetrize in the four variables. The connected functions 
On, are then obtained by (|4.3|) and ([4.9|), and they are simply given by: 



On{\u . . . An) = J^^;^ {—) ^n(Al+zO,...,An + zO). 

Universality. We observe that the only universal features of Ci(A) = p(A) are the 
square-root singularities at the edge of the distribution, otherwise the function is potential 
dependent. Instead the two-point correlation function has a universal form, and can thus 
be calculated from the gaussian model. 

The n-point smoothed correlation functions with n > 3 can be calculated recursively 
by a systematic method described above. The main property is that they consist in a sum 
of a finite number of universal functions and involve only the first n — 2 moments of the 
potential. That means for instance, that two potentials V and V* induce the same three- 
point function as soon as they yield the same Mi and Ji coefficients, but they don't need 
to be identical. However, the n-point correlation function is no longer given by a gaussian 
model, since it is of order l/A^^"^"^. Perturbative corrections to the gaussian model have 
to be considered. 

Note that the determination of correlation functions allows an evaluation of the mo- 
ments Mfc, in a case where the potential V is unknown. Let us emphasize that these 
moments play an important role in the study of critical points. Since 

p(A) = ^M(A)V(A-a)(6-A) 

we see that if some moments vanish, the behaviour of p near the end-points is no longer a 
square-root. When the m first M^ vanish, one finds p ~ (A — a)"^"'"^/^ which corresponds 
to the w}^ multicritical point of the one-matrix model of 2D gravity. [^^. 

12 



5. Conclusions 

In this article we have recovered, by a completely different method, the results of 
ref. |jl[] concerning the two-point eigenvalue correlation function of a random hermitian 
matrix in the limit in which the size A^ of the matrix becomes large. Our method has 
allowed us to generalize the results to new distributions. In addition, we have established 
a recursion relation between successive n-point correlation functions at leading order for 
A^ large. 

Brezin and Zee in ref. [IHhave shown that the two-point function is universal, and 
therefore identical to the function of the gaussian matrix model. The gaussian model 
is the fixed point of a renormalization group p|J^] and a direct RG analysis should be 
performed to put this result in perspective. 

In the same way, Brezin and Zee have shown that higher correlation functions vanish 
at leading order. The contributions we have calculated here should be considered as 
corrections to the leading scaling behaviour. The explicit expressions we have obtained 
show that they depend now on successive moments of the potential, indicating an implicit 
classification of the deviations from the gaussian model in terms of their irrelevance for A^ 
large. Here also it would be interesting to confirm the qualitative aspects of these results 
by a direct RG analysis. 

Appendix 1. Orthogonal polynomials P^i An explicit expression 

Let us show that the orthogonal polynomials P^ defined by the orthogonality condi- 
tion: 



are given, up to a normalization, by equation 



,2 -mrV(M) 
Pr,iX)oc (TMe ^ ^det(A-M). 



First, this integral clearly yields a polynomial of degree n in A. Let us then verify the 
orthogonality property: after integration over the unitary group, the integration measure 
J dM reduces to an integration over the eigenvalues of M, and the Jacobian of this trans- 
formation is a square Vandermonde determinant: 

PniX)ocf n dA,e~^^^'^'\A-AO A'(Ai,...,An), 

1=1, ...,n 

13 



where 

A(Ai,...A„) = n(Az-A,). 



i<j 



In this form we recognize a more classical expression [^ . Then, setting A = Aq : 

<Pn(Ao)-A^>oc /" n dA,e~^^^^'^A(Ao,Ai,...,A„)A(Ai,...,A„)AS^. 

i=0,...,7i 

The first Vandermonde is completely antisymmetric in the n + 1 variables, we can therefore 
antisymmetrize the factor A(Ai, . . . , An)Ao', the result is zero if m < n because the only 
polynomial completely antisymmetric and of degree less than n — 1 in Aq is zero. Thus 

< p,^ . A"^ >= . 

it proves that P^ is orthogonal to any P^ with m < n, and that < Pn-Pm >= as soon 
as 771 7^ n. 

Remark that a similar integral representation can also be found for multi-matrix 
models. Consider two families of orthogonal polynomials Pn and Qn, such that: 

f f -N{V{X) + U{fi)-cXfi) ^ ^ 

<Pn{X)-QmW>= / dAd/Xe PnWQmW OC 6n7n- 

Then, we have: 

f f -NtvlViMi) +U(M2) - CM1M2] 

P^(A)= / ZdMidMae L v 1; v 2; ^^ det(A - MJ 



where Mi and M2 are hermitian n x n. Similarly: 

f f -Ntv[V(Mi) + U(M2) - CM1M2] 

Q„,{^)= ZdMidMse ^ ^ ^' ^ '' '^det(/i-M2) 



Appendix 2. Connected correlation functions and the kernel k(A, fi) 



Some exact expressions exist for the correlation functions of matrix models fll 
terms of a kernel k(A, h). We again consider the matrix distribution ( |1 . 1| ) : 

, , 1 -NtvV(M) 
V{M) = -e 

14 



m 



The corresponding measure can be rewritten in terms of the eigenvalues A^ of M and a 
unitary transformation U which diagonahzes M: 

V{M)dM = Z-^dU Yl dA,e~ ^ '^ A'^{Xi, . . .,Xn) 

i=l...N 

(A = ni<7 (-^j ~ ^j) being the Vandermonde determinant). Therefore, the probabihty that 
the eigenvalues of M are Ai, . . . , Xn is: 

PAr(Ai,...,AAr) oc JJ e ' A'^ {Xi, . . . , Xn), 

i=l...N 

and the correlation functions are obtained by partially integrating over some eigenvalues: 

p(Ai) = / Y\ dXiPN{Xi,...,XN), 

•^ i=2,...,N 

P2(Ai,A2)=/ JJ dXiPN{Xi,...,XN) 
^ i=3,...,N 

... and so on. 

The Vandermonde determinant A can be written as: 

A(AO = detAr\ 
and thus, after some linear combinations of columns of the matrix: 

A = det n,_i(AO, n^(A) = x^ + o (A"-i) , 

identity true for any set of polynomials {n„} normalized as above. In order to perform the 
A integrations, we choose 11^ oc P„, P„ being the orthogonal polynomials (|2.1|) . A^ is the 
product of two such determinants, therefore it is the determinant of a matrix product: 

A^ ocdeti^(Ai,Aj) (2.1a) 

N-l 

K{X,, A,) = J2 Pk{X^)Pk{XJ). (2.16) 

fc=0 

The proportionality constant in (|2.1aD is here irrelevant because the eigenvalue distribution 



is normalized. The Darboux-Christoffel formula (appendix 3) tells that: 

Pn{X)Pn-M - Pn{p)Pn-i{X) 



K{X, n) = a- 



A — // 
15 



where a is a normalization constant depending on A^, and a = {a — b)/4 when A^ -^ oo. 
The important properties of K{X, y) are: 

due ^^K{X,u) = l, due ^ ^ K{u,u) = N , 

-NV(u) 
due ^ ' K{\,u)K{u,y)=K{\y). 

An exphcit expression for p^ 

1 /■ TT TT -NV{\i) 

Pn(Ai,..., An) = Try / _[_[ dAi YY e del K{\i,\ 



J7 7 

i=n+l,...,Ar i=l,...,Ar 



can be obtained by successively integrating over eigenvalues. Using the rules of K integra- 
tion it is easy to prove by induction 



A^"(A^-n)! 



Pn (Al, A2, . . . , \n) = ^, det K (Ai, Xj) , 



where we have introduced the reduced function 



1 -{N/2)[V{X) + V{p)] 
K{X,p) = —e K{X,p). 



Therefore, when the A,- are all distinct we have: 



»z 



(C»(Ai) . . . 0(An)) = det k(A„ a,) = ^(-1)'^ n ^(-^^' -^-^) 

The connected function will involve only the sum over the permutations a such that 
n^ ^(Ai, A(jJ cannot be split in the product of two cyclic products, i.e. only cyclic per- 
mutations will contribute to the connected function. This intuitive result is a classical 
combinatorial identity, can be proven fromThe connected function can thus be written: 



permutations a i 

- (-1)^+1 [k (Al, A2) /. (A2, A3) . . . « (An, Al) + ■ ■ ■] 



where the additional terms in the r.h.s. symmetrize the expression. Or in a more compact 

way: 

n 

(o(Ai) . . . o(An))_ = i-ir+' Yl n ^ (^- ^-^) 

cyclic permutations a i=l 
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Appendix 3. Derivation of the Darboux— ChristofFel formula 

The polynomial XPn{X) can be expanded on the basis of the Pm with m < n + 1: 

n+l 

XPn{X) = 2_^ QnmPm\X)- 
m=^n — l 

The orthogonality condition ( p.l|) implies that the matrix Q is symmetric: 
The polynomial (A — n)K{X, /i) can thus be written: 

Af-l / N N 

fc=0 \i=0 j=0 

All the terms cancel, except the upper-bounds: 

(A - fi)K{X, y) = Qn,n-i {Pn{X)Pn-M - PN(;u)Piv-i(A)), 

therefore a = Qn,n-i- In the large A^ limit, it is possible to calculate a. The simplest 
way of doing this is to calculate APjv-i from expression ( p.6|) . Since 

cos If 1 COS {^jJ — (f) = cos {ip — if) — (cos ^jJ + cos {^jJ — 2if)) , 



we have: 

,^ a + b ^ b — a , ^ ^ 

XPn-1 = -^-Pn-1 - -^- {Pn + Pn-2) 

and therefore a = (a — b)/A. 
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